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Gravitational Waves from a Fissioning White Hole 
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We present a fully nonlinear calculation of the waveform of the gravitational radiation emitted in 
the fission of a vacuum white hole. At early times, the waveforms agree with close approximation 
perturbative calculations but they reveal dramatic time and angular dependence in the nonlinear 
regime. The results pave the way for a subsequent computation of the radiation emitted after a 
binary black hole merger. 
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INTRODUCTION 



The omputation of gravitational radiation from the in- 
spiral and merger of binary black holes poses a difficult 
boundary value problem. In the geometrically simplest 
and physically most natural treatment, the black holes 
are modeled by the gravitational collapse of a pair of stars 
(or other astrophysical bodies). However, this is a chal- 
lenging hydrodynamic problem which requires simulating 
a pair of orbiting bodies for a sufficient time to verify a 
negligible amount of incoming radiation in the initial con- 
ditions, then following their subsequent collapse to black 
holes and finally computing the outgoing radiation in the 
exterior spacetime. Alternatively, in the purely vacuum 
approach, the individual black holes form from implod- 
ing gravitational waves. This avoids hydrodynamical dif- 
ficulties at the expense of a globally complicated initial 
value problem. The imploding waves may emanate either 

(i) from a past singularity or (ii) from past null infinity 
X~. In case (i), the appropriate boundary condition at 
T~ is that there be no ingoing radiation but, assuming 
the time reversed version of cosmic censorship, the past 
singularity implies a white hole horizon *Hr on which 
boundary data must be specified in some arbitrary man- 
ner in order to determine the exterior spacetime. In case 

(ii) , ingoing radiation from X~ is present at early times 
when the black holes are formed but ingoing radiation 
must be absent at late times in order for the outgoing 
radiation to be unambiguously attributed to the merging 
black holes. In this work, we present a solution to the first 
stage of a new two-stage global treatment of the vacuum 
binary black hole problem g). The approach, based 
upon characteristic evolution, has been carried out in the 
regime of Schwarzschild perturbations where advanced 
and retarded solutions of the linearized problem can be 
rigorously identified || . Computational experiments are 
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necessary to study the applicability of the approach to 
the nonlinear regime. 

From a time- reversed viewpoint, this first stage is 
equivalent to the determination of the outgoing radiation 
emitted from the fission of a white hole in the absence 
of ingoing radiation. This provides the physically correct 
"retarded" waveform for a white hole fission, were such 
events to occur in the universe. Although there is no 
standard astrophysical mechanism for producing white 
holes from a nonsingular matter distribution, white holes 
of primordial or quantum gravitational origin cannot be 
ruled out. 

This fission problem has a simpler formulation as a 
characteristic initial value problem than the black hole 
merger problem. The boundary of the (conformally com- 
pactified) exterior spacetime contains two null hypersur- 
faces where boundary conditions must be satisfied: past 
null infinity X~ , where the incoming radiation must van- 
ish, and the white hole event horizon Tl~ , which must de- 
scribe a white hole, which is initially in equilibrium with 
no ingoing radiation and then distorts and ultimately fis- 
sions into two white holes with the emission of outgoing 
gravitational waves. If we approximate T~ by an outgo- 
ing null hypersurface J~ , which intersects Tl~ at an early 
time (approximating past time infinity i~) close to the 
initial equilibrium of the white hole, then data on these 
two null hypersurfaces, TL~ and J~ , constitute a stan- 
dard double-null initial value problem, whose evolution 
determines a portion of the exterior spacetime extend- 
ing to T + , where the outgoing radiation is computed. In 
contrast, the corresponding problem for the "retarded" 
waveform from a black hole merger involves two disjoint 
null hypersurfaces where boundary conditions must be 
satisfied: past null infinity X - , where the incoming ra- 
diation must vanish, and the future event horizon 7i + , 
which describes the merger of the two black holes and 
their subsequent approach to equilibrium. 

In previous work j2j, we treated the fission problem 
in the close approximation pj as a perturbation of a 
Schwarzschild background. In this paper we present a 
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fully nonlinear treatment that reveals new and interest- 
ing strong field behavior. We carry out the evolution of 
this vacuum double-null problem by means of a charac- 
teristic evolution code j|[ [| [?J , using a recent version of 
the code which improves accuracy in the highly nonlinear 
region || . Caustics in the ingoing null hypersurfaces used 
to foliate the exterior spacetime restrict the evolution to 
the pre-fission stage. 

We use a conformal horizon model (|, [lC], [ll) to sup- 
ply the necessary null data for a horizon corresponding 
to a white hole fission. The conformal horizon model 
provides a stand-alone description of the intrinsic null 
geometry of the horizon. The algorithm for generating 
horizon data is constructed to handle a general event 
horizon representing the fission of a spinning white hole 
into two outspiraling white holes of non-equal mass [ p"T| . 
The specific application in this paper is to the axisym- 
metric head-on fission into equal mass white holes. (The 
necessary data and evolution codes are, however, not lim- 
ited to the axial symmetry of a head-on collision.) The 
resulting horizon geometry is an upside-down version of 
the standard trousers-shaped event horizon for a binary 
black hole merger in the time-reversed scenario. 

We study a range of models extending from the per- 
turbative close limit, in which the fission occurs in the 
infinite future, to the highly nonlinear regime. Nontriv- 
ial global changes, accompanied by dramatic time depen- 
dence of the horizon geometry, arise in passing from the 
perturbative to the highly nonlinear regime 11. The ex- 
istence of a marginally trapped surface divides the hori- 
zon into interior and exterior regions, analogous to the 
division of the Schwarzschild horizon by the r — 2M bi- 
furcation sphere. In passing from the perturbative to the 
strongly nonlinear regime there is a transition in which 
the fission occurs in the portion of the horizon visible 
from I + . Thus these results reveal two classes of bi- 
nary white hole spacetimes, depending upon whether the 
crotch in the trousers (where the fission occurs) is bare, 
in the sense that it is visible from T + , or hidden inside a 
marginally trapped surface. In this paper we evolve this 
data using the characteristic code to study the properties 
of the gravitational radiation produced by this dramatic 
behavior of a white hole fission. 

In Sees. || and |ll we review the formalism and data 
necessary for a characteristic evolution of the spacetime 
exterior to a dynamic white hole by means of the char- 
acteristic code. In Sec. IV we present the detailed wave- 



forms for a one-parameter family of spacetimes varying 
from the close approximation to the highly nonlinear 
regime in which the fission is visible from X + . 

We retain the conventions of our previous papers 
^, 0, |^, |l0) , with only minor changes where noted in the 
text. For brevity, we use the notation f iX — d x f to denote 
partial derivatives and / = d u f to denote retarded time 
derivatives. We represent tensor fields on the sphere as 
spin-weighted variables jl2| in terms of a dyad qA for the 
unit sphere metric qAB — <Z(a<7b) m some standard choice 



ical code uses two overlapping stereographic coordinate 
patches.) We compute angular derivatives of tensor fields 
in terms of S and 9 operators |ll|, e.g. f t A = 5R(?a9/), to 
compute the gradient of a spin- weight zero (scalar) field / 
in terms of the spin-weight 1 field 3/ and the spin- weight 
-1 field §/. 



II. THE DOUBLE-NULL PROBLEM FOR A 
FISSIONING WHITE HOLE 

We treat the fission of a white hole by a double null 
initial value problem based upon the white hole horizon 
and on an outgoing null hypersurface 3~ , which em- 
anates from an early time slice S~ of H.~ approximating 
the initial equilibrium of the white hole. The horizon 
pinches off in the future where its generators either caus- 
tic or cross, producing the (upside-down) trousers picture 
of a fissioning white hole. 

The double-null problem, first formulated by 
Sachs [jnj, is most conveniently described in Sachs 
coordinates consisting of (i) an affine null parameter 
u along the generators of H~, which foliates H~ into 
cross sections S u and labels the corresponding outgoing 
null hypersurfaces J u emanating from the foliation, (ii) 
angular coordinates x A which are constant both along 
the generators of T~t~ and along the outgoing null rays 
of J u , and (iii) an affine parameter A along the outgoing 
rays normalized by V"uV a A = — 1, with A = on H~ . 
In these x a = (u,X,x A ) Sachs coordinates, the metric 
takes the form 



ds 2 



— (W - g AB W A W B )du 2 - 2dudX 
2g AB W B dudx A + g AB dx A dx B (2.1) 



of spherical coordinates, e.g. x 1 



(The numer- 



We set qab = r 2 fiAB, where det(/i J 4s) = det(gAs) = 
q(x A ), with qAB the unit sphere metric. We represent 
the conformal metric \iab by the complex spin- weight 2 
field J = \h,ABq A q B ' ■ The remaining dyad component, 
given by the real function K = ^hABQ A Q E \ is fixed by 
the determinant condition K 2 = 1 + J J. 

The requirement that the horizon be null implies that 
W = on H~ . In addition, we fix the gauge freedom 
corresponding to the shift on H" so that n a d a = d u is 
tangent to the generators, implying that W A — on H~, 
The choice of lapse that u is an affine parameter implies 
d\W = on 7i~ , We further fix the affine freedom by 
specifying u = u- on the early slice S~ approximating 
the asymptotic equilibrium of the white hole in the past. 
The outgoing null hypersurface J~ emanating from S~ 
approximates past null infinity X~ . 

The affine tangent to the generators of Hr , n a d a = d u 
satisfies the geodesic equation n b Vfcrt a = and the hyper- 
surface orthogonality condition n' a V fc rt c l = 0. Following 
the approach of Refs. [l(|, we project four-dimensional 
tensor fields into the tangent space of TC~ using the op- 
erator 
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P" 



(2.2) 



where l a = — V a it. 

The evolution proceeds along the outgoing null hyper- 
surfaces J u emanating from the foliation of TL" and ex- 
tending to (compactified) X + . In this problem, the com- 
plete (and unconstrained) characteristic data on TL" are 
its (degenerate) intrinsic conformal metric Hab{u, x a ), 
or cquivalcntly J(u,x A ), expressed in terms of the afhnc 
parameter u. In addition, the characteristic data on J~ 
are its intrinsic conformal metric /iab(A, x a ), or J(A, x A ) 
expressed in terms of its affine parameter A. 

The remaining data necessary to evolve the exterior 
spacetime consist of the intrinsic metric and extrinsic 
curvature of S~ (subject to consistency with the char- 
acteristic data) jl4|, [19] . This additional data consist of 
the surface area r, the inward expansion r, the outward 
expansion r t \ and the twist [fl6| 



of initial conditions on S" . The value of r on TL" de- 
termines the convergence of the ingoing null rays. Once 
the intrinsic geometry J and the area coordinate r are 
known, the vacuum equation Rau — is used to prop- 
agate the twist u) — q A LOA along the generators of TL~ . 
This determines lo in terms of its initial value at S~. 
The Rab = vacuum equations propagate the outward 
expansion and shear of the foliation S u along TL" . The 
trace part propagates the outgoing expansion determined 
by r t \. The trace-free part is an evolution equation for 
h-AB,\, which describes the shear of the outgoing rays. 

In summary, the data for the double null problem in- 
clude the conformal metric Iiab on TL" and J" and the 
quantities r, r, w and r t \ on S" . Equations (Al), (2.29), 
(2.30), and (2.34) of Ref. [jll| propagate the data on S" to 
all of TL~ . The data required on S~ can be inferred from 
the asymptotic properties of the white hole equilibrium 
at i~ . The propagation equations, given in spin- weighted 
form in [ p"l"| , are implemented numerically with a second 
order Runge-Kutta scheme jl8| . Complete details can be 
found in Ref. 111. 



u a = P*n c VJ b = (0,0, ua) 



(2.3) 



(Our choice of shift implies that uj u = lo\ = 0.) The twist 
is an invariantly defined extrinsic curvature property of 
the u = const cross sections of TL" , independent of the 
boost freedom in the extensions of n a and l a subject to 
the normalization n a l a = — 1. 

We use the conformal horizon model to supply data 
J on TL~ corresponding to a fissioning white hole for a 
sequence of models ranging from the clo se a pproximation 
to the highly nonlinear regime (see Sec. III). On J" , we 
set J = to model the absence of ingoing radiation. In 
carrying out the evolution computationally, the first step 
is topropagate the data on S~ along the generators of 
TL ]l7] so that it can be supplied as boundary data for 
the exterior characteristic evolution code. 



B. Bondi-Sachs coordinates and characteristic code 
variables 

The null code is based upon the Bondi-Sachs version 
of the characteristic initial value problem |l9|, ^o) . It is 
designed to evolve forward in time along a foliation of 
spacetime by outgoing null hypersurfaces. The Bondi- 



Sachs coordinates differ from the Sachs coordinates (2.1 



by the use of a surface area coordinate r along the out- 
going cones rather than the affine parameter A. Because 
a generic horizon is not a hypersurface of constant r, it is 
advantageous to first discuss the necessary data in terms 
of Sachs coordinates and then transform to the r coor- 
dinate. In Bondi-Sachs variables, the metric takes the 
form 



A. Propagation equations on the horizon 

Einstein's equations for the double-null problem de- 
compose into (i) hypersurface equations intrinsic to the 
null hypersurfaces J u , which determine auxiliary met- 
ric quantities in terms of the conformal metric Hab] (ii) 
evolution equations which determine the rate of change 
duhAB of the conformal metric of J u \ and (iii) propaga- 
tion equations which are constraints that need only be 
satisfied on TL" H]. The Bianchi identities ensure that 
the propagation equations will be satisfied in the exterior 
spacetime as a result of the hypersurface and evolution 
equations. Integration of the propagation equations de- 
termines the additional horizon data necessary for the 
characteristic evolution. 

One of the propagation equations is the ingoing Ray- 
chaudhuri equation R uu — 0, which propagates the sur- 
face area variable r along the generators of TL" in terms 



ds 2 = - (e 2 ^ - r 2 h AB U A U B ) du 2 - 2e 2 Hudr 

~2r 2 h AB U B dudx A + r 2 h AB dx A dx B . (2.4) 

The field U A is represented in spin-weighted form U = 
q A U A . 

The evolution of the exterior spacetime requires a 
transformation of the data on TL" from Sachs coordi- 
nates (u, A, x a ) to Bondi-Sachs coordinates (u,r,x a ), as 
described in detail in Ref. (ll). Since r is not in general 
constant on TL~ , the horizon does not lie precisely on ra- 
dial grid points r^. Consequently, an accurate prescrip- 
tion of boundary values on the n grid points nearest the 
horizon requires a Taylor expansion of the horizon data. 
Restricted to TL" , the metric variables r and J have the 
same values in both Sachs and Bondi-Sachs coordinates 
and our choices of lapse and shift imply that the Bondi- 
Sachs variables (3, U and V are related to Sachs variables 
by Eqs. (2.15), (2.31), and (2.32) of Ref. 
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A similar construction jTI| provides their first r deriva- 
tives in terms of known Sachs variables. We obtain 
J r from Eq. (2.33) of Ref. fn|, /3, r from the R\\ = 
Raychaudhuri equation for the outgoing null geodesies, 
Eq. (2.37) of Ref. @, and the value of d r U on the hori- 
zon from the auxiliary field Q = q A QA, where 



r 2 e 



!0 h AB U 



(2.5) 



Here Q is obtained from the twist and other Sachs vari- 
ables as per Eq. (2.39) of Finally, we compute d r V 
by obtaining V \ from the A derivative of Eq. (2.32) of 
Ref. @. 

The values of each metric function (J, j3, U, V) and its 
first radial derivative can be used to consistently initialize 
field values at the 7Vgrid points near the horizon. Bound- 
ary values for the code are then provided, to second order 
accuracy, on grid points bracketing the horizon. For each 
outgoing null ray the value of the area coordinate r-j-i- 
on the horizon is known. This value is bracketed by the 
nearest grid points rj_i < < r^, where the metric 
values are computed by the second order accurate Taylor 
expansion 



J r , = Ah- + (n ~ t h-)JAh- 



(2.6) 



After the metric quantities have been computed at the 
grid points neighboring the horizon, the code can evolve 
the spacetime exterior to the horizon. 



C. Evolution equations 

The system of equations that determines the exterior 
spacetime forms a hierarchy pi] , |22j . Alternative formu- 
lations are available; see for instance Ref. Q. Here we 
use a recent formulation [^J, which was specifically de- 
veloped to handle the extremely nonlinear post-merger 
regime of binary black hole collisions. It reduces all an- 
gular derivatives to first order by introducing the aux- 
iliary variables v = 9J, k — <5K and B = 9/3. This 
results in the hierarchy of hypersurface equations and 
one complex evolution equation for the conformal met- 
ric function J given by Eqs. (16)-(25) of Ref. ||. We 
note that the terms on the right hand side of the evo- 
lution equation, Eq. (16) of have been grouped into 
hypersurface terms Jjj which vanish for linear perturba- 
tions of a spherically symmetric spacetime and a term 
P u which isolates the only nonlinear term containing a 
(retarded) time derivative of J. The introduction of the 
auxiliary variables B, v, and k in |^] eliminates all second 
angular derivatives from the hypersurface and evolution 
equations. This leads to substantially improved numeri- 
cal behavior plover the standard characteristic formula- 
tion |^|, ||, p3| |24| in the extremely nonlinear post-merger 
regime of binary black hole collisions. Details are given 
in Ref. B. 



We retain the radial and time integration schemes of 
the standard characteristic formulation. The radial in- 
tegration algorithm is explained in detail in Refs. || M, 
B^ , ^djj . In a departure from Ref. 0L we use a three-step 
iterative Crank-Nicholson scheme ]27[ to ensure stability 
of the time evolution Pq| . 

The evolution code first integrates the hypersurface 
equations [Eqs. (17)-(23) of J§|]] on the initial null hy- 
persurface J~ in the region exterior to the horizon. The 
evolution equation, Eq. (16) of ||] is then used to ad- 
vance to the next hypersurface in the region exterior to 
the horizon, etc., advancing the computation of the ex- 
terior spacetime in the region from the horizon to X + . 
The evolution continues as long as the coordinate system 
remains well behaved. 



D. The Bondi news 

The Bondi news function N(u,x B ) jL9|, [|| |(| is an 
invariantly defined field on X + which gives the amplitude 
of the flux of gravitational wave energy. Its computation 
requires evaluation of the conformal factor Q con f neces- 
sary to compactify the spacetime in a conformal Bondi 
frame |n| , which is an asymptotic inertial frame in which 
the slices of I + have unit sphere geometry. Following 
the approach of Ref. we set fl con f = w con //r, where 
(u, x ) is a smooth non- vanishing field at J+. At 
the early time w_ at which we initiate the evolution, 
the conformal metric Hab of the outgoing null hyper- 
surface approaches the unit sphere metric qAB at X + so 
that LOconf- = w co „/| M _ = 1. However, as the evolution 
proceeds, the (u, x A ) coordinates, which are naturally 
adapted to H~, become non-inertial at I + . This results 
in the time dependence 



{d u + L A d A )\og 



D A L A 



(2.7) 



where L A is the asymptotic value of U A at T + || . 

In addition, a transformation to conformal Bondi co- 
ordinates (uB,y A ) on X + is necessary to determine the 
Bondi news 



M(u B ,y A ) = N(u(u B ,y A ),x A (uB,y A )), 



(2.8) 



as measured in an inertial frame. The relevant expres- 
sions are given in Sec. IV and Appendix B of Ref. In- 
spection of Eqs. (B1)-(B6) of Ref. § reveals that second 
angular derivatives of the conformal factor u) con f enter 
the calculation. For improved accuracy, we remove these 
second derivatives by introducing the auxiliary vari- 
able W = QtUconf- Since w con / is defined solely on X + , 
we use the consistency relation >V „ = Buj CO nf,u to prop- 
agate W along the generators of T + , initializing it by 
W_ = 5w con /_. We use a combination of second-order 
Runge-Kutta and mid-point rule schemes || for the time 
integration of Eq. (2.7). 
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III. NULL DATA FOR THE AXISYMMETRIC 
HEAD-ON FISSION 

The conformal horizon model §, || supplies the con- 
formal metric Hab constituting the null data for a binary 
black or white hole. For the case of a head-on fission of 
a white hole, the model is based upon the flat space null 
hypersurface Ti emanating from a prolate spheroid Sq of 
eccentricity e and semimajor axis a, which is embedded 
at a constant inertial time t = in Minkowski space. 
Traced back into the past, Ti expands to an asymptoti- 
cally spherical shape. Traced into the future, Ti pinches 
off at points where its null rays cross or at caustic points 
where neighboring null rays focus. 

With the dyad choice q A — (l,£/sin#), the intrinsic 
conformal metric Iiab of Ti as a flat space null hypersur- 
face is determined by the spin-weight-2 field 



J 



q A q a h AB {t 1 x°) 



















r<p/ 






re J 



(3.1) 



where rg and are the principle radii of curvature of 
the spheroid. The metric of the spheroid induced by its 
embedding in a flat space is c/ab — f 2 hAB, where r 2 sin 9 
is the local surface area. 

The white hole horizon shares the same manifold 
Ti and the same (degenerate) conformal metric as its 
Minkowski space counterpart and Eq. ( |3.1[ ) provides the 
conformal null data for the white hole horizon in the t 
foliation. But the surface area and affine parametrization 
of the white hole horizon and its flat space counterpart 
differ. As a white hole horizon, Ti extends infinitely far 
to the past of Sq to an asymptotic equilibrium with fi- 
nite surface area. The intrinsic metric of the white hole 
horizon is given by qab — Q 2 gAB, and has local surface 
area r 2 sin# where r = fir. The conformal factor fl is 
designed to stop the expansion of the white hole in the 
past so that the surface area radius asymptotically hovers 
at a fixed value r — Rqq. The conformal factor is given 

by@0 



which governs the second derivative of the surface area 
with respect to the affine parameter. Because the focus- 
ing power determined by their conformal geometries are 
the same, the Raychaudhuri equation implies 



U 2 r = Idlr. 



(3.3) 



Consequently, the different behavior of the surface area 
coordinates r and f due to the conformal factor f2 implies 
that the affine parametrization t on the white hole hori- 
zon is related to its flat space counterpart u according 

to m 



dt 
di 



(12u(u-p) - a 2 ) 2 

(5 P + m-2m) 2(2p/ai+1) 
(5p-^-2tt) 2 ( 2 P/A'- 1 ) 



where 



v/l3p 2 - a 2 



(3.4) 



(3.5) 



and where the affine scale is fixed by the condition 
dt/dt — > 1 as t — > — co. The scale of the particular affine 
parameter u used for the time step in the evolution code 
is related to t by u = t(l + ecos 2 9) 3 ^ 2 (chosen so that 
the null data J has early time behavior of the pure spin- 
weight 2 qua drup ole form J ~ — esin 2 9/u). 

Equation (3.4) determines the rate of deviation of a 
slicing adapted to an affine parameter u of the white hole 
horizon from the t slicing given by the Minkowski embed- 
ding. The build up of a large relative angular dependence 
between the u and t foliations leads to the change in 
topology of the fissioning white hole and the associated 
pair-of-pants shaped horizon. We begin our evolution at 
an early slice S~, at u = u_, where the horizon is close 
to equilibrium as a Schwarzschild white hole. Initially, 
the u and t foliations behave quite similarly. The time 
at which nonlinear behavior becomes significant is con- 
trolled by the location of the Minkowski spheroid Sq, at 
t = 0, with respect to the startup slice S~ . 



where p is a model parameter, a = tq — t^ is the difference 
between the principal curvature radii, and u is a flat space 
affine parameter along the generators of Ti with the same 
scale as t but with its origin shifted along each ray by 
u = i — h( r 6 + r <j>)- For an initially Schwarzschild white 
hole of mass M, — 2M. Smoothness of the white 
hole requires that the parameter p > om/V^3, where 
<jm is the maximum value of a attained on Sq. (For the 
prolate spheroid considered here, <jm = ae.) 

Both as a flat space null hypersurface and as a white 
hole horizon, Ti must obey the Raychaudhuri equation, 



IV. WAVEFORMS FROM FISSIONING WHITE 
HOLES 

We use the characteristic code [|| to evolve the exterior 
spacetime of an axisymmetric, head-on white hole fission 
and compute the detailed waveforms for a sequence of 
spacetimes in the range < e < 10~ 2 , varying from the 
perturbative to the highly nonlinear regime. The object 
is to search for new gravitational wave physics resulting 
from the highly nonlinear behavior of a binary horizon. 
For that purpose, the simulations reported here were car- 
ried out on a grid of 77 x 77 x 137 (stereographic patch 
x radial) points. At this resolution, a higher eccentricity 
evolution (e ~ 10 -2 ) requires approximately one hour on 
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a single Pentium II 750 Mhz processor. Since the nu- 
merical code has been demonstrated to be second order 
convergent, more accurate results, especially concerning 
late time behavior near the point of fission can be at- 
tained by increasing the numerical grid resolution. Even 
though the characteristic code [|| is not parallel, param- 
eter search runs of the type described here are entirely 
feasible on computing clusters with the present code. 

We initialize the simulations according to the specifi- 
cations in Ref. where the horizon geometry for this 
data was computed and studied. We begin the evolution 
at u = u_ = —100 with the initial mass M_ = 100 and 
conformal model parameters i_ = — 10 (which fixes the 
location of the spheroid So), a — 1 (the semi- major axis 
of the spheroid) and p = (10~ 2 + 10~ 5 )/\/T3 (just above 
the minimum value of p allowed by regularity of the con- 
formal model for this range of eccentricities jll| . On J~ 
(u = m_), we set J = to model the absence of ingoing 
radiation in an initially Schwarzschild white hole. 

The limit e — ► yields the Schwarzschild solution and 
small eccentricity corresponds to the close approximation 
where the fission takes place far in the future behind a 
marginally trapped surface so that it is hidden from ex- 
ternal observers. For sufficiently small e, the close ap- 
proximation is valid in the entire exterior Kruskal quad- 
rant but for large e it is valid only at very early times. 
For large e, the fission occurs in the region of spacetime 
visible from I + . For a critical value e — e c 10~ 4 there 
is a transition between these two regimes. For e < e c , 
the evolution terminates when the expansion of the out- 
going null hypersurfaces vanishes along some ray. (In the 
Schwarzschild limit this occurs simultaneously on all rays 
on the black hole horizon.) For e > e c , the evolution ex- 
tends to the point of fission when r — ► on the equator 
of the white hole. 

Figure |l| shows the time dependence of the fractional 
mass loss AM(u)/M_ = [M(u) - M_]/M_ for a se- 
quence of runs. A negligible amount of mass is radiated 
at startup, confirming that in all cases the initial geom- 
etry is very close to Scharzschild. The mass loss remains 
small until u w —90, near the location of S$ on 7i~, at 
which time a rapid change in dt/du occurs, and the mass 
loss rises sharply. 

Following the sharp rise in AM/M_ , a transient period 
is observed in all the curves. The largest contribution 
to the radiated energy comes from the last stage of the 
evolution. The total mass loss remains well below the 1% 
level until eccentricities of about e 10~ 5 are reached, 
when the total mass loss at the end of the simulation 
approaches 2% of the initial white hole mass. 

At higher eccentricities, i.e. at e > 10~ 4 , most of the 
initial mass is radiated away during the later stages of 
the simulation. The rate of mass loss in the final stage is 
particularly dramatic in the cases e > e c where the fission 
is visible from X + . 

Figure |2| shows again the mass loss profiles of Fig. [I], 
but now scaled by the relative amplitude (e/eo) 2 , with 
€q = 10~ 6 . The agreement seen clearly indicates that the 
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FIG. 1: Relative mass loss AM /Mo as a function of the hori- 
zon affine parameter u, shown here for values of the eccen- 
tricity parameter ranging from e = 10 -6 to e = 10 -2 . 
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FIG. 2: The relative mass loss AM/M , scaled by e 2 and 
overlayed on the profile for the case e = 1CT 6 . The striking 
agreement indicates that the mass loss behavior is still 0(e 2 ). 



behavior is 0(e 2 ) in the range where they overlap. 

Figure |^ shows the time dependence of the news mea- 
sured at the equator for e = 10~ 6 to e = 10 -2 . Here we 
have also overlayed the plots, rescaling them by (e/eo), 
with eo = 10~ 6 . We plot the real part of the Bondi 
news at a point on the equator, as a function of Bondi 
time u b ■ Due to axisymmetry, the imaginary part of the 
Bondi news is zero to within discretization error. The 
sharp pulse clearly visible in the news at ub ~ — 52 re- 
sults from the choice of parameters in the white hole 
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FIG. 3: The real part of the Bondi news at a point on the 
equator, scaled by e and overlayed on the profile for the case 
e = 10 -6 , as a function of Bondi time ub- A sharp pulse is 
clearly visible at ub ~ — 52, which is a distinct feature of the 
white hole horizon data. 



FIG. 4: The real part of the Bondi news at early times 
(u — —96.63) for e = 10 -6 . At this early time, the angu- 
lar dependence is purely quadrupolar. 



horizon data which control the rapidity of the fission. 
Here, in order that the model yield a bare fission, the pa- 
rameters have been chosen so that the fission occurs on 
an extremely rapid time scale. Subsequent to this sharp 
pulse, the news undergoes a damped oscillation (which 
is too early in the final ringdown to be associated with a 
quadrupole quasinormal mode). 

The Bondi news as computed by the full code osculates 
the results obtained in the perturbative limit |B2] , in the 
regime in which the perturbative and nonlinear codes are 
clearly solving the same problem (up to ub ~ 100). At 
later times (ub ~ 105), the fully nonlinear news calcu- 
lation deviates from the perturbative results as a con- 
sequence of the appearance of higher harmonics. These 
higher harmonics, which are clearly visible in the news as 
a function on the sphere, arise from the nonlinearity of 
the equations and cannot be observed in a perturbative 
evolution. 

Figures [| and || illustrate the angular behavior of 
the news at early time (u — —96.63) and at late time 
(u = —68.64), respectively. We plot the real part of the 
Bondi news on a single stereographic patch for the case 
e = 10~ 6 . The graphs clearly show that at early times the 
news is pure quadrupole, in agreement with the perturba- 
tive regime of the close approximation. The nonlinearity 
of the problem subsequently introduces higher harmonics 
in the news. These higher harmonics are clearly visible 
at later times in the plots of the Bondi news as a function 
on the sphere. The most notable feature is the dimple 
seen in Fig. [| on the profile of the news at the pole (the 
center of the stereographic patch) . 

This feature is partially due to the retardation effect in- 



RelNcws] 




FIG. 5: The real part of the Bondi news at late times (u — 
—68.64) for e = 10~ 6 . The dimple visible in the angular profile 
is due to the presence of higher harmonics, I > 2. 



troduced by the angular-dependent redshift of the Bondi 
frames at future null infinity X + . We have chosen the 
initial cut of X + , corresponding to a constant affine time 
u on the horizon, to be a cut of constant Bondi time. 
We then follow the inertial observers at X + to define a 
Bondi time slicing. Figure which displays the angular 
dependence of Bondi time ub for a fixed horizon time u, 
reveals a pronounced increase in Bondi time at the pole 
relative to the equator. 
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FIG. 6: The Bondi time as a function on a stereographic patch 
at null infinity at a late times (it = —59.92) for e = 10 -6 . The 
pole is "running ahead" of the equator, leading to the dimple 
observed in Fig. H. 



sented here are for the axisymmetric non-spinning case, 
the data and evolution codes are not restricted to any 
symmetry. It will be interesting to see how the results 
for a head-on collision are modified in the fission of a 
spinning white hole. 

Reexpressed in terms of the time-reversed scenario of a 
black hole merger, the boundary conditions for a fission 
corresponds to no outgoing radiation in the black hole 
case. Nevertheless, the results pave the way for an appli- 
cation of characteristic codes to calculate the fully non- 
linear waveform emitted in a binary black hole collision 
in the time period from merger to ringdown. Waveforms 
from a black hole merger can be expected to differ from 
those from a white hole fission, as has been observed in 
close approximation studies The fission process is 
directly observable at X + whereas the merger waveform 
arises indirectly from the preceding collapse of the matter 
or gravitational wave energy that forms the black holes. 
This suggests that the fission is a more efficient source 
of gravitational waves and that the high fractional mass 
losses computed here cannot be attained in a black hole 
merger. 
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